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The work described in this report was performed by the 
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Abstract 

Definite integrals are evaluated for the cardinal 
functions of an interpolation method which provides 
continuity over a triangular grid. 
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1. Introduction 


In Ref. 1, algorithms are given foi C^-corapatible in- 
terpolation over a plane triangle given values of f, f > and f 

X / 

at the vertices of the triangle. In particular on p. 30 of 

Ref. 1, functions and y. are defined such that with 

B- -8 /2 and y.=y./2 the interpolated value w can be ex- 
11 11 

pressed as 


Cl) w := E (f.a.+f -B'+f -Y-) 

1 1 x.x^i y,i 


In this note we derive expressions for the definite 
integrals of the cardinal functions a^,B^, and i*l,2,3 
over the triangle. These integrals are needed, for example, 
if a surface representation problem includes conditions on 
volumes to be encompassed under the surface. 


2 . Preliminary integration formulas 


From Ref. 2, we obtain the following formulas for in- 
tegration of polynomials over a triangle: 


(2) -CPx “ A Pl^i’3’3^ 

(3) -CP 2 “ ^ 

(4) /P 3 = ^^C25p2[-5-,j,j] + 25p2fj,^,^] + 25p2[-g-,j,'^]-27pj(j,2,^]} 


1 
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In the above formulas the integration is over a triangle in the 
xy plane. The integrand denotes a polynomial of degree at 
most i in X and y. The symbol A denotes the area of the triangle. 
The notation PiCrj,r 2 »rj] denotes the value of the polynomial 
Pj^ at the point in the triangle whose barycentric coordinates 
are r^^, T 2 > and r^. For discussion of the affine invariance of 
quadrature formulas stated in terms of barycentric coordinates 
see Ref. 2. 

The above formulas can be used to obtain the following 
integrals for the fundamental linear, quadratic, and cubic func- 
tions used in Ref. 1. 


(5) 

^"i 

" A/3 

(6) 

/<l>i 

= A/12 

(7) 




In the following two sections we determine for 

the two different definitions of p ^ given in Ref. 1. Formulas 
for the integrals of the cardinal functions are then derived 
in Section 5. 

3. Determining for the rational function 

In Section 7, pp. 22-23, of Ref. 1 p^ is defined 
as the Zienkiewicz rational function 

‘’i * i-1,2,3 

Consider the particular triangle Tq having vertices = (0,0), 
^2 ® (1>0), Vj = (0,1). In this triangle we can relate car- 
tesian xy coordinates to barycentric coordinates as follows: 

r^ = 1-x-y 

t2 = X 


r3 = y 


2 
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Then 


. 2 . .2 


= X y (1-x-y)/ t(l-x) (l-y)l 


Expand the term (1-x-y) as (l-x)+ (1-y) -i 
sum of the following +hrcc integrals: 


Then / pn 

T 

^0 




0 


0 


lo = 


^ y^dy dx 

0 ^ 0 


I, * - 


“ -|S 


2 2 
1 X^ ,1-x yL. 


1-x i 


1-y 


dy dx 


Evaluating I^: 




I |f^(y^-2y^*/) dy 


Ifl . 2 + 3, _ 1_ 
4 5^ 90 


The integral I 2 has the same value as I^. The inner 
of Ij is 

^ dy = * 2(l-y) - £n(l-y)}^‘'' 

= + 2x-£nx + j- 2 

» -j(x^-4x+3) - tn X 

» -j(x-l) (x-3) - In X 



is the 


integral 


3 
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Thus 


^3 “ "I ^^x^(x-3) dx + dx 

= 1 + fl ^ 

8 4 1-x 

_ 3 ^ y.1 In X - /^ £n X dx - x£nx dx 

8 0 1-^ 0 0 

3 TT^ r /> ^ » x^ 1 

= j ^ - [xtnx-x + j- Inx - ^]q 

2 

= I - J- + I » (39-47t2)/24 

Combining results we obtain: 

( «) i »1 ■ II * h * '3 - W * 5!f ^ " (593-60f2)/360 

^0 

Since the area of triangle Tq is 1/2, the mean value of pj^ 
over Tq, and thus over any arbitrary triangle, is (539-60v^)/180. 
Furthermore, since the integrals of p^^,P2» and pj are the same 
we may drop the subscripts and simply write the mean value of 
any of these three functions as 


( 9 ) P = (593-60ti^)/180 

I 0.00457 63107 47992 


The author thanks Dr. E.W. Ng for doing the analytic in- 
tegrations of this section. The result given in Eq. ( 9 ) has 
been corroborated by a numerical integration of Pj^ over the 
triangle Tq which gave agreement in the first seven significant 
digits. See Section 8 for more details on ths numerical integration. 


4 




4. Determining /p^ for the piecewise cubic function 

In Section 8, pp. 24-26, of Ref. 1 , is defined 

for i=l as the Clough-Tocher piecewise cubic function 


"i 


f riiarjrj.rjCSrj 


3)1/6 if 

if 
if 


rj^ * min{rj^,r2»r j} 
r2 * minirj^,r2,r5} 
rj = min{rj^,r2,rj} 


Thus p^ is defined as a different cubic polynomial in each of 
the three subtriangles indicated in Fig. 1. 




The common intersection point of the three subtriangles has bary- 
centric coordinates • 

Within each subtriangle the integration formula, Eq. ( 
can be applied. The evaluation points for this formula in sub- 
triangles Tj^ and T 2 are indicated in Fig. 2. 


V 
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1 

Table 1. Integration of over [ 6 r 2 Tj+rj^ (Srj^-3) ] 


Point 

Barycentric 

Coordinates 

Value of p^ 

Multiplier 
[See Eq. ( 4) ] 

Product 

a 

1 2 2 
5, 5, 5 

7 

3-125 

25 

189/(81-5) 

b 

1 10 4 

15, 15, 15 

4 

81.5 

25 

100/(81*5) 

c 

1 4 10 

15, 15, 15 

4 

81-5 

25 

100/(81*5) 

d 

14 4 

9, 9, 9 

37 

81 -27 

1 

-27 

-185/(81*5) 

— 

Sum: 

204/(81 *5) = 63/(27'*5) 


1 2 

Table. 2. Integration of Pj^ over T 2 . Pj^ = ^ "’^2'*’^^3^ 


Point 

Barycentric 

Coordinates 

Value of p^ 

Multiplier 
[See Eq. ( 4)1 

Product 

e 

2 12 
5, 5, 5 

1 

2 *3 *25 

25 

135/(2 -81 *S) 

f 

4 1 10 

29 

25 

29/(2 -81 -5) 

15, TT, 15 

2 *81 *125 

g 

10 1 4 

1?, TI, I5 

11 

2 *81 *175 

25 

11/(2 *81 *5) 

D 

4 14 

9, 9, ^ 

11 

2 *3^ 

-27 

-55/(2 *81 *5) 


Sum: 

120/(2*81*5)= ^ 


6 
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f 

1 


The arithmetic involved in applying Eq. (4) to in sub- 

triangles and T 2 is summarized in Tables 1 and 2. The 
integral over T^ will be the same as that over T 2 * 


Let A denote the area o£ the whole triangle, 
each subtriangle T^ has area A/3. It follows that 


/Pi - 


A 


68 

27-5 


2.4 

27 


180 


Then 


The integrals /P 2 and /p^ have the same value as fp^. Thu: 
the mean value of any one of these functions p^ is given by 


( 10 ) 


p = 1/180 = 0.0055555555. . . 


The result given in Eq. (10) has been corroborated by a num- 
erical integration of Pj^ which gave agreement in the first 
seven significant digits. See Section 8 for more details on the 
numerical integration. 

5. Integrals of the cardinal functions and 


1 :; 


Referring to p. 30 of Ref. 1, introduce B = jB 

1 - ^ 
and Y = jY so that the final equation for w can be written 

simply as 




\ 


Use a bar over the symbols p^, q^, a^, B^^. and 
to denote the average value of each of these functions over the 
triangle. For example, ^ -^Pi’ case of the functions 

pj the integral over 
we simply write p- * 
we have ^ = (593-60 tt 
function and p •» 1/180 if pj^ is the Clough and Tocher piecewise 
cubic function. 


the triangle is independent of i. Thus 
^/p^ for i»l,2,3. From Eq. (9)and(10) respectively 
2)/180 if p^ is the Zienkiewicz rational 


7 





33-808 


Let 

\ !^\ 

Then, using the notation and equations of Ref. 1, p. 30, and 
Eqs. (5-7) of the present paper, we obtain 


/8i 

It 

3X^/p^ = 3a. 

Pi 

= 

^^iP l2 

^i 

= 

3X-P - 1 

^ 12 

^i 


I ^ 3(X._, 

^i 


‘ I 

^i 

= 

I ^^-lPi-1 ^ 


Finally we have 
(11) /w . A 

3 

Note, as a matter of interest, that E a. = 1, independent 

— i= 1 ^ 

of the value of p . 


6 . Example; An isosceles right triangl e 

Consider the isosceles right triangle of Fig. 3. 

whose essential geometric parameters are given in Table 3. This 

2 

triangle has area A - a /2. 



Figure 3 


(a,0) 
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Table 3. Parameters of the triangle of Fig. 3. 


f' 

i 

^i 

^i 

h 

'i 

r 

f 

1 

-a 

a 


0 














2 

0 

-a 

a 

-1 


3 

a 

0 

a 

1 


Choosing the functions to be the Clough-Tocher 

piecewise cubic functions we have p - 1/180. Using the form- 
ulas of Section 5, we obtain the values in Table 4. 


Table 4. Mean values of auxiliary and cardinal functions 


i 

Pi 

‘li 

“i 


^i 

1 

5/60 

-5/60 

22/60 

6a/120 

6a/120 

2 

4/60 

-6/60 

19/60 

-9a/120 

5a/120 

3 

6/60 

-4/60 

19/60 

5a/120 

-9a/120 


These coefficients should provide exact quadrature for 

any quadratic function. As a specific example let a*120 and 

2 

consider the function f(x,y) » x . 

By the formula of Eq.( 3 ) the integral of f over 
the triangle in question is given by 

(12) I - I {60^ + 60^} » 10-120^ 

9 
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Alternatively, usinii Table 4, we note that only two 

2 

of the nodal values are non- zero, namely f (120,0) = 120 and 
f^(120,0) = 240. Thus using Eq.(ll) and 0^2 and 82 from 
Table 4, we compute the integral 

T 120^ [120^-19 ^ 

^ = — [—60 - 5 - 240 J 

= 10 - 120 ^ 

which is in agreement with lq.(12). 


7. Example; An equilateral triangle 


Consider the equilateral triangle of Fig. 4, whose 

essential geometric parameters are given in Table 5. This 

1/2 2 

triangle has area A = 3 a /4. 



Figure 4 


Table 5. Parameters of the triangle of Fig. 4. 


i 

“i 




1 

-a/2 

3^'^^a/2 

1 

0 

2 

-a/2 

-3^/^a/2 

1 

0 

3 

a 

0 

1 

0 


10 
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Since all X. are zero the functions will not 
enter into the problem. Using the formulas of Section S, we 
obtain the values in Table 6. 

Table 6. Mean values of auxiliary and cardinal functions 


i 

Pi 

^i 

“i 

^i 

^i 

i 

1/12 

- 1/12 

1/3 

a/16 

3^/^a/48 

2 

1/12 

- 1/12 

1/3 

-a/16 

3^/^a/48 

3 , 

1/12 

-1/32 

1/3 

0 

j 

-2 -3^^^a/48 


These coefficients should provide exact quadrature for 

any quadratic function. As a specific example let a = 16 and 

2 

consider the function f(x,y) = y . 

By the formula of Eq.(3 ) the integral of f over 
the triangle in question is given by 

(13) I = j {(4-3^/^)^ + (4-3^/^)^} 

= 32A = 2^^ • 3^/^ 

Alternatively, using Table 6, we note that only two of 

1/2 

the nodal values are non-zero, namely f(8,8’3 ' ) = 3*64 and 
fy(8,8 -3^/^) = 16 •3^'^^. 

Thus using Eq.(ll) and oTj and Y3 from Table 6, we 
compute the integral 


11 
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T . A r3-64 16-3^/^-2-3^/^*16, 

J ' A {-r-3 ^ } 


* 32A 


= 2 II . .1/2 


which is in agreement with Eq.(13)* 


t 

1 
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8. Corroboration of results by numerical integration 

As a check against possible blunders in the determination 
of values of p in Sections 3 and 4, direct numerical integration 
Was used to obtain estimates of p for the two different definitions 
of p. 


The quadrature subroutine used was the JPL Univac 1108 
library subroutine RMB1/RMB2 (Ref. 3) . This subroutine effects 
quadrature over a two-dimensional region (a triangle in the present 
case) by making an appropriate sequence of calls to the library's 
one -dimensional quadrature subroutine ROMBS/ROM2. This system of 
subroutines was developed by W. R. Bunton and M. R. Diethelm of 
JPL using an adaptive algorithm based on Romberg quadrature. 

Results of these numerical integrations are summarized 
in Tables 7 and 8. For each of the two functions being studied, 
the analytically determined value of p is known. For each of the 
two function^, the numerical integration was done five times with 
five different values of the accuracy tolerance parameter. 

Letting p denote the value obtained from the i^ numerical inte- 
gration, values of the difference p^^ - p are given in Tables 7 § 8. 

The results of these numerical integrations are consistent 
with the analytically determined values of p. 


13 



33-808 


Table 7. 


Numerical integration of the Zienkiewicz rational 
function. 

p = (593-60 it^)/180 * 0.00457 63107 47992 


i 

Absolute 

accuracy 

requested 

Number of 
evaluations 
of integrand 

- P 

1 

10'^ 

113 

-0.00000 14834 

2 

10"* 

113 

-0.00000 14834 

3 

io-=’ 

215 

+0.00000 01164 

4 

10’^ 

263 

+0.00000 00059 

5 

10"^ 

1413 

+0.00000 00007 


Table 8. Numerical integration of the Clough-Tocher piecewise 
cubic function. 

p » 1/180 = 0.00555 55555 ... 


1 

Absolute 

accuracy 

requested 

Number of 
evaluations 
of integrand 

P^ - P 

1 

XO- 

133 

-0.00000 31391 

2 

10'* 

141 

-0.00000 59111 

3 

10'^ 

231 

-0.00000 00555 

4 

10"^ 

681 

+0.00000 01064 

5 

lo'^ 

895 

-0.00000 00055 


14 
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